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I. INTRODUCTION 

Noise-induced escape from a metastable state underlies 
many important phenomena in nature and technology. It 
is often investigated experimentally for systems close to 
a bifurcation point where a metastable state disappears. 
Examples include Josephson junctions 0,0,11], nanomag- 
nets 000], and mechanical nano- and microresonators 
0013- Studying escape in this regime is advantageous 
in two respects. First, near a bifurcation point the ac- 
tivation barrier is small, so that the escape rate can be 
measured in a conceivable time even for low temperature 
or small noise intensity. Second, the dynamics near a 
bifurcation point displays system-independent features, 
including the onset of one or several soft modes. In turn, 
this leads to scaling of the escape rate with the control 
parameter of the system. 

We will discuss escape for system near a generic saddle- 
node bifurcation, which is one of the most common types 
of bifurcations. Here, for very small values of the con- 
trol parameter r\ counted off from the bifurcation value 
7]b = 0, the metastable and unstable states of the sys- 
tem are close to each other in phase space. The motion 
along the axis that connects them is slow. This is easy 
to understand when metastable and unstable states are 
a minimum and a saddle point of the potential energy 
of a particle. The curvature of the potential has oppo- 
site signs at these points. Therefore the potential locally 
flattens out as rj approaches zero. In agreement with 
this picture, for very small rj the system dynamics is con- 
trolled by a single soft mode q, which is overdamped. 

In this paper we consider systems with very small 
damping. Their motion may become underdamped even 
where the parameter r\ is still close, although not too 
close to the bifurcation value. The dynamics of the soft 
mode is characterized by the vibrational frequency at the 
metastable state oj a . It is much less than other frequen- 
cies, and lu u — > for T) — > 0. Yet u a may exceed the 



relaxation rate, which is characterized by an indepen- 
dent small parameter e. Such a situation often occurs in 
underdamped Josephson junctions, for example. 

Noise-induced escape near a bifurcation point is well 
understood in the limit 77 — s- 0, where the motion is over- 
damped mill. Here, under fairly general as- 
sumptions the activation energy of escape scales as if'/ 2 . 
This agrees with the scaling of the mean-field free-energy 
barrier near a termination point of the metastable state 
[l4| . As we show, such scaling holds also for under- 
damped systems further away from 77 = provided the 
metastable and saddle states are still close to each other. 
This can be understood for a Brownian particle in a po- 
tential well. Here, the activation energy of escape is given 
by the height of the potential barrier, independent of fric- 
tion; friction only affects the prefactor in the escape rate 
[l5j . Therefore the scaling is the same in the under- and 
overdamped regimes. The scaling rff 2 was discussed al- 
ready in the early work on Josephson junctions [Lsl ] and 
later for underdamped driven oscillators |17| and mag- 
nets m- 

As we show, besides the if/ 2 scaling, underdamped 
systems may display a different scaling near, but not too 
close to the bifurcation point. This is because the mo- 
tion slowing down is not necessarily related to the stable 
and unstable state approaching each other. The slow- 
ing down may have a different, nonlocal origin. It can 
be thought of as resulting from flattening, or an overall 
decrease in magnitude of the effective potential of the 
system in a broad range of the generalized coordinate. 
We provide a theory of escape for this nonlocal case and 
find the corresponding critical exponent. The results are 
applied to a specific system, which provides an example 
of this scaling behavior and has attracted much attention 
recently 001113. 

In Sec. II we discuss the model of a fluctuating under- 
damped system. In Sec. Ill we describe two scenarios 
that lead to the motion slowing down and ultimately to 
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the saddle-node bifurcation, in the presence of friction, 
as the control parameter 77 — > 0. The effect of friction 
is described in Sec. IV. Energy diffusion is considered in 
Sec. V. In Sec. VI we obtain the critical exponents for 
the activation energy of escape, which is the central re- 
sult of the paper. In Sec. VII the theory is applied to 
a nonlinear oscillator driven by a strong resonant field 
and compared with the previous results for this model 
[r^.[T^| . Sec. VIII contains concluding remarks. 



II. THE MODEL 

Close to a saddle-node bifurcation point one of the 
motions in the system becomes slow compared to other 
motions. Quite generally it is described by a Langevin 
equation 

q = d p H{p,q)-ev^(p,q)+f^(t), 

p = -d q H(p, q) - ev^ip, q) + / W (t). (1) 

Here, q and p are the coordinate and momentum of the 
slow motion, and H is its effective Hamiltonian. The 
terms ev^ q ' p ) and f( p ' q ^ describe friction and noise, re- 
spectively. We have assumed that the motion is slow on 
the scale of the correlation time of the reservoir coupling 
to which leads to friction and noise. Therefore we dis- 
regard retardation in the friction force. In the same ap- 
proximation we disregard the noise correlation time and 
set the noise to be (5-correlated in time, 

</«(*)/«(*')> =2£>«(p, «)<*(* -iO, (2) 

where i,j enumerate the q and p components of noise. 
The diffusion matrix is symmetric, D qp = D pq , and gen- 
erally Dij depends on the dynamical variables q,p. The 
noise intensity oc maxDy will be the smallest parameter 
of the theory. 

We assume that the friction coefficient e is small and 
that the friction force is non-Hamiltonian. The latter 
requires that 

Vv = d q v M + d p v^ (3) 

is not equal to zero. The part of the force (v^ q \v^) 
for which Vv = can be incorporated into the Hamil- 
tonian, leading to its small renormalization cx e. This 
concerns, in particular, ^-independent terms in v^ q > and 
p-independent terms in v( p >. 

The Hamiltonian H (p, q) = H(p, g; 77) depends on the 
control parameter 77 which characterizes the distance to 
the bifurcation point. This point is located at 77 = 0, for 
e = 0. We require that (i) H{p, q; 77) be analytic in 77 in a 
neighborhood of 77 = 0, and (ii) d v H ^ for 77 = 0. These 
two conditions determine the order of the Hamiltonian 
H with respect to 77 uniquely for the different functional 
forms of H used below, i.e., if both rji and 772 meet them, 
then 772 = 0(7yi). 



The physical picture of motion leading to escape is 
simple. The system is initially prepared at some point 
(q,p) in the basin of attraction to the metastable attrac- 
tor (focus) (q a ,p a )- Over the relaxation time t T cx e , 
the system approaches the attractor. Noise leads to small 
fluctuations about the attractor, but for small noise in- 
tensity the system stays close to it for a time that largely 
exceeds t T . Ultimately there happens a large fluctuation 
that carries the system over the basin boundary leading 
to escape from the metastable state. 

A. Brownian particle near a saddle-node 
bifurcation point: standard analysis 

The probability W of noise-induced escape from a 
metastable state is a complicated function of the dynami- 
cal parameters. However, it displays a universal behavior 
in the vicinity of a bifurcation point. This behavior has 
been well understood in two limiting cases. The best 
known case is a Brownian particle in a potential well, 
which escapes from the well due to thermal fluctuations. 
The motion of a Brownian particle of mass m = 1 in a 
potential U{q) is described by the equation with 

H(j>,q) = ^p 2 + U(q), fM=0, 

v (i) = 0, wW = p, (4) 

and with f^ p \t) being white noise of intensity D = 2ekT. 

The dynamics of a Brownian particle is special in sev- 
eral respects. First, the value of the momentum in a 
stationary state is independent of the parameters and is 
equal to zero. Second, independent of the friction coeffi- 
cient, the stable and unstable states of the system q a and 
q s are the local minimum and maximum of U(q). Third, 
as a consequence of the form of the noise, the escape rate 
has a simple form 

W = const x e - AC//feT , AU = U{q s ) - U{q a ). (5) 

The exponential factor here is just the ratio of the Boltz- 
mann factors in the populations of states with energies 
of the saddle point and the attractor. It is independent 
of the friction coefficient. In contrast, the prefactor in 
W depends on e [l5j : however, in this work we will be 
interested only in the exponent in W. 

The potential U(q) depends on the control parameter 
77. The standard analysis of this dependence goes as fol- 
lows. When 77 is close to the critical value 77 = 0, the 
stationary states are close to each other and ultimately 
merge for 77 = 0; we set q a = q s = for 77 = 0. For small 
77 the potential is a cubic parabola, U(q) « t4ut>(<z), 

U cuh {q) = -q 3 /3 + 7 N (6) 

locally in q (which has been appropriately scaled). The 
stable and unstable states are q a = —i] 1 ^ 2 , q s = 77 1 / 2 . 
Both the distance q s — q a and the frequency of vibrations 
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about the potential minimum in the absence of friction 
uj a = (477) 1 / 4 decrease with decreasing 77, as expected. 

Eq. JBJ) leads to the scaling of the activation barrier 
with 77 of the form AU = (4/3)?/ 3 / 2 . We emphasize that 
this scaling holds independent of friction, for a Brownian 
particle. Because the friction coefficient e is small, the 
approximation © applies both where the system is un- 
dcrdamped, with u a 3> e, or undcrdamped, with tu a <C e. 
The critical value of the bifurcation parameter 77 is also 
independent of e. 

B. Close vicinity of the bifurcation point 

Irrespective of the friction coefficient, very close to the 
bifurcation point the dynamics of the system is over- 
damped [2£j. It is described by a soft mode Q, which is 
the distance between the saddle and the attractor. The 
Langevin equation for Q after proper scaling and chang- 
ing to a slow time r can be written as 

§ = Q 2 -v + f (Q) (r), (7) 

For 77 > the system has a stable and unstable states 
Q = =R/^ 2 , which merge together for rj = 0. 

Because the motion near the bifurcation point is slow, 
with characteristic time ry -1 ^ 2 , and occurs in a narrow 
range \Q\ < 77 1 / 2 , the noise f^'{r) can be considered 
white, with a Q-independent intensity D. Then the es- 
cape rate is [T(| 

,,1/2 A 

W= 1 ^e- R ' D 1 R=-r?l 2 . (8) 

7T 3 

Eqs. 0, JSJ can be understood as the limit of Brown- 
ian motion for s 3> 77 1 / 4 , where the Brownian motion 
becomes overdamped. Eq. (JJJ follows from Eqs. 0, (@J 
in this limit, with Q — q and t = et. The activation 
energy of escape R © scales with 77 in the same way as 
the barrier height AU for the corresponding Brownian 
particle. 

We emphasize that the applicability of Eq. (J7J) is not 
at all limited to systems close to thermal equilibrium. 
It applies essentially to any weakly fluctuating dynami- 
cal system close to a saddle- node bifurcation point, see 
[l3j and references therein. The scaling R oc if 1 / 2 has 
been considered a benchmark of activated escape near 
a saddle-node bifurcation point. It has been used very 
broadly, in particular as a tool for measuring critical cur- 
rent of Josephson junctions 

III. SADDLE-NODE BIFURCATION IN 
UNDERDAMPED NONEQUILIBRIUM SYSTEMS 

We will consider the general case of a system away from 
thermal equilibrium, which is described by Eq. with 
small friction coefficient e. Our major assumption is that, 



as we decrease the control parameter of the Hamiltonian 
77 from a positive value ~ 0(1) while keeping e =const, 
the system will experience a saddle-node bifurcation. Be- 
cause the system is nonequilibrium, the bifurcation value 
r\B may depend on the friction parameter e. For conve- 
nience, in what follows we do not set t\b = 0, but we 
assume 

< 1, lim 775 = 0. 

For small s the system is underdamped for t/b <C 
?7 <C 1. Still very close to r/s it invariably becomes 
overdamped and is described by Eq. Q. We want to 
find the forms of the Hamiltonian and the correspond- 
ing phase portraits of the Hamiltonian trajectories that 
will be compatible with this scenario of the saddle-node 
bifurcation. 

By assumption, the bifurcation results from the change 
of the Hamiltonian dynamics, not from the special struc- 
ture of the friction force. Then a necessary condition is 
that the Hamiltonian system has a center (q c ,Pc) and a 
saddle point (q s ,p s ) for 77 > 0. With weak friction, the 
center becomes a focus. When 77 becomes very close to 
rjB the focus becomes a node and ultimately merges with 
the saddle point. 

Central to the occurrence of a bifurcation is the slowing 
down of the system with decreasing rj. We note that this 
slowing down should occur not only near the center, but 
in the whole range of energies between the energies of 
the center E c = H{q c ,p c ;rj) and the saddle point E s = 
H(q s ,p s ;r]). We will discuss two possible phase portraits 
that meet this condition. They are shown in Fig. ^ 



A. Local Hamiltonian bifurcation 

The simplest evolution of the phase portrait of Hamil- 
tonian dynamics corresponds to the center and the sad- 
dle point approaching each other in phase space with 
decreasing rj and ultimately merging for 77 = 0, sec upper 
panel in Fig. ^ We call this "local Hamiltonian bifur- 
cation" . The frequencies of the Hamiltonian trajecto- 
ries surrounding the center uj{E) decrease with 77 — > 0, 
because these trajectories are pressed against the homo- 
clinic trajectory which has zero frequency. 

This case is easy to analyze. We expand the Hamilto- 
nian H in q,p counted off from their value at the degen- 
erate point into which the center and the saddle point 
merge for 77 = 0. Because of the degeneracy, the determi- 
nant d 2 Hd 2 H — (d q d p H) 2 = at this point. We choose 
q along the zero-eigenvalue eigenvector of the matrix of 
second derivatives of H for 77 = 0. Then to leading order 
in q,p, rj the expansion of H takes the form 

H(p,q;v) « \p 2 + U cnh {q) (9) 
with U cuh (q) = U cuh (q;rj) given by Eq. ©. 
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FIG. 1: Two types of phase portraits leading to a saddle-node 
bifurcation in the presence of dissipation. The homoclinic tra- 
jectories go through the saddle point (q a ,p s = 0). The center 
(q c ,Ps — 0) [attractor, in the presence of dissipation] is shown 
by the cross. As the control parameter of the Hamiltonian r\ 
approaches the critical value, the loop in the left panel shrinks 
in all directions. The center and the saddle approach each 
other in phase space and ultimately merge at the bifurcation 
point. The frequency of vibrations about the center u a de- 
creases, because the center becomes closer to the homoclinic 
trajectory, which has an infinite period. In the right panel, 
the homoclinic orbit also shrinks as r\ — > 0, but the stationary 
states remain separated. Still the frequency ui a decreases with 
decreasing 77, because the trajectories are pressed against the 
homoclinic orbit. 



In Eq. wc appropriately scaled q,p and disregarded 
all cubic terms except for the terms oc q 3 and all terms 



oc rj except for r/q. An estimate \q\ 



3/4 



1/2 



\P\ 



u{E) 



oc rj 



1 ' 4 shows that the higher-order terms 



in q,p,ri lead to small corrections for 77 <C 1. In or- 
der to make the results more intuitive, we chose the 
coefficient at p 2 to be positive. In this case the cen- 
ter q c = — r/ 1 / 2 , p c — has lower energy than the saddle 
point q s — 77 1 ' 2 ,p s = 0. If we chose H = — p 2 /2 + U cu b(q) 
it would be the other way around, as we illustrate on the 
example below, but this does not change the following 
arguments. 

In the approximation @ the homoclinic loop that sur- 
rounds the center has a simple form 



p = ± 



1/2 



(10) 



The vibration frequency about the center is to a = (477) 1 / 4 . 
As expected, uo a — > for t] — ► 0. The difference of the 
energy values in the extreme points 



B. Nonlocal Hamiltonian bifurcation 

Slowing down in the Hamiltonian system occurs also 
if the homoclinic trajectory shrinks down with rj — ► 
while the center (q c ,p c ) and the saddle point (q s ,p s ) do 
not move or move only slightly, sec lower panel in Fig.^ 
As in the local case, the frequencies of all trajectories 
inside the homoclinic loop go to zero for 77 — > 0. We call 
it a nonlocal Hamiltonian bifurcation. It corresponds 
to the homoclinic orbit becoming degenerate, with no 
motion along it at all, at any point. The values of the 
Hamiltonian at the center and the saddle point coincide, 
H(q c ,p c ;V = 0) = H{q s ,p s ;?i = 0). 

We choose the momentum p as the variable along 
which the trajectories are shrinking for 77 — >• 0. This vari- 
able is limited by the "height" of the homoclinic loop in 
Fig-ffl an d therefore it is small for small Since on the 
degenerate homoclinic trajectory q = d p H = for 77 = 0, 
after a proper rescaling of variables the Hamiltonian for 
small 77, p takes the form 



H{p,q;rj) = ^p 2 +rjU{q). 



(12) 



A more general form 



H = 



a?(«) 



p-i] 1/2 a 2 {q) 



V U(q) 



is reduced to Eq. (|12|) by a transition to canonical vari- 
ables p' = 01(g) [p— 7] 1/2 a 2 (q)],q' = J dq/ai(q) (we as- 
sume that ai does not have zeros in the range of interest). 

The specific form of the potential U (q) is not impor- 
tant, except that U (q) must have a local minimum and 
maximum. For d 2 H > [as chosen in Eq. (TSjl]. the 
minimum of U(q) corresponds to the center q c , and the 
maximum corresponds to the saddle point q s . The homo- 
clinic trajectory goes around q c , and therefore the turn- 
ing point of this trajectory q t (where p = 0) lies on the 
opposite side from q c with respect to q s , as in Fig. ^ 
We assume, in agreement with the qualitative picture of 
approaching the bifurcation that, in the interval [qt,q s ], 
the potential has no other extrema besides q c and has no 
singularities for real q . 

A simple form of U(q) consistent with the above con- 
ditions is a cubic parabola, U(q) = — q 3 /3 + q. In this 
case q c = — l,q s = 1, and the turning point q t = — 2 
independent of ?y. However, since the distance between 
g c and q s does not decrease with decreasing 77, there are 
no reasons to believe that the potential is gencrically a 
cubic parabola. We will give an example of a different 
potential below. 

For the Hamiltonian (|12|) , both the momentum and fre- 
quency scale with 77 as \p\ oc uj(E) oc tj 1 ^ 2 . As expected, 
Lj a = r^^JJ' '(q c ) — > for 77 — > 0. To leading order in rj 
the energy difference at the saddle point and the center 
is 



AE = H(q s ,p s ;r,) - H(q c ,p c ; V ) = -ij 3 ' 2 . 



(11) 



AE = H(q s ,p s ;r))-H(q c ,p c ;ri) 

= V [U(q s )-U{q c )]. (13) 
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It follows from Eqs. (|13fl . that for the local and 

nonlocal Hamiltonian bifurcations AE displays scaling 
behavior with the control parameter, AE oc rfi. However, 
the scaling exponents £ are different, £ = 3/2 and £ = 1 
in the local and nonlocal cases, respectively. 



IV. EFFECT OF FRICTION 

The analysis of friction is particularly simple in the 
case of a local Hamiltonian bifurcation 0. Generally, 
one would expect the friction force v^ p,qS> to be smooth 
near the point q = p = 0. For small \p\ it can be 
expanded to linear terms in q,p. Keeping only the terms 
that are not reduced to renormalization of the parameters 
of the Hamiltonian, we have 



The center (q c 



a„q, 



^ 1//2 jPc = 0) becomes a stable state 



,(p) 



a p p. 



(14) 



provided e(a q +a p ) > 0. It shifts because of the friction to 
(q a = —(?? + £ 2 a p a q / A) 1 / 2 , p a = ea q q a ). The bifurcation 
value of the control parameter in the local case is 



T7 loc 
Vb 



—£ 2 a p a q /4:. 



(15) 



The analysis of the nonlocal case is somewhat more 
involved. In this case friction changes the behavior of the 
system very significantly: while the center and the saddle 
point remain separated for rj — > in the Hamiltonian 
case, because of the friction the corresponding attractor 
and the saddle point are approaching each other as rj goes 
to the bifurcation value rf^. 

To understand the constraints on the friction force we 
note that far from the bifurcation point, rj 3> rj^j, the 
motion of the system is slowly decaying oscillations. It is 
characterized by the period-average value Ve of the rate 
of energy change 



E = -e(v {p) d p H 



v^O q H) 



It is convenient to do the averaging by writing the integral 
over time as an integral over the trajectory H (jp, q; rj) = 
E, which gives 



V E 



ME) 
u(E) 



dqv (p) - dpv {q) 

dqdpVv, (16) 

A(E) 



where the divergence of v is defined by Eq. ©. The 
circulation is taken along the trajectory H (p, q; rj) = E, 
and A(E) is the enclosed region. 

We assume that the surface integral in Eq. (|l(i|> does 
not change sign for all energies between the center and 
the saddle point, E c and E s . If the integral is posi- 
tive (we assume e > 0) then the period-average energy 
monotonously drifts down to its value E c if it was ini- 
tially between E c and E s . This is consistent with the 



assumption that friction does not lead to extra station- 
ary states. A sufficient and physically plausible condition 
is that Vv > C > everywhere in the area enclosed by 
the homoclinic trajectory in Fig. ^ with C independent 
of e,r]. We note that Vv must be positive at the point 
(qa,Pa) in order for the state (q a ,Pa) to be stable. 

Close to a nonlocal Hamiltonian bifurcation point 



\P\ 



1/2 



<C 1, and therefore, gencrically, we can ex- 



pand v(p, q) in p. The term oc p in can be disre- 
garded, and it suffices to keep only a linear in p term in 
i/p); the p-independent term in can be incorporated 
into the Hamiltonian leading to a trivial renormalization 
U{q) -> U{q) + (E/r))v^(p = 0,q). Then the stationary 
states in the presence of friction are given by the equa- 
tions 



p = ev^\ 
vd q U{q) + e 2 



(17) 



v {q) d p v ip) 



p=0 



The positions of the states remain practically unchanged 
with the varying parameter of the Hamiltonian rj as long 
as rj ^ e 2 . When this condition holds, the frequency of 
vibrations near the attractor uj a ~ rj 1 / 2 largely exceeds 
their decay rate oc e and the motion is undcrdamped, as 
expected. 

If, as we assume, the system displays a saddle-node 
bifurcation for a small rj, then with decreasing rj the sta- 
tionary states come closer to each other and ultimately 
merge together for rj = ifg ~ e 2 . Operationally, the 
bifurcation value of the control parameter r/g is deter- 
mined by the condition that two solutions of Eq. l(T7|) 
for q,p coincide. An illustration for a specific model is 
provided below. 

We note that Eq. (|17f) may have two solutions for rjg. 
This means that there are two bifurcation points in the 
region of small \tj\. The possibility of such behavior is 
clear from the Hamiltonian l|12|) . For example, for U(q) of 
the type of a cubic parabola, the patterns of trajectories 
for positive and negative rj are similar, with the center 
and the hyperbolic point interchanged. The two values 
of r;^ 1 correspond to the saddle- node bifurcations of these 
two types of motion with decreasing \rj\ in the presence 
of friction. 

Eq. i|16[) describes energy drift in the case of a lo- 
cal bifurcation as well. In this case, from Eq. i|14f> . 
Vv = a.p + a q is constant. The integral (|16f> can be 
found explicitly, in terms of elliptic functions. 



V. ENERGY DIFFUSION 

A major effect of noise on the dynamics of under- 
damped systems is drift and diffusion of energy. Their 
rates are determined by the noise intensity parameters 
Dij J2J). In the weak-noise limit, where Dij are the small- 
est parameters of the theory, noise-induced energy drift 
can be disregarded compared to the drift induced by fric- 
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tion. Therefore we will be interested only in energy dif- 
fusion. 

The period- average value De of the energy diffusion 
coefficient can be obtained from the equation of motion 
in a standard way: we write dH/dt as d q Hq + d p Hp, 
substitute the noise terms and for q and p, 



find the correlator of dH/dt using the noise correlators 
Pjl. and average the corresponding diffusion coefficient 
over the period 2tt/lu(E). Changing at the last step 
from time integration to integration along the trajectory 
H(p, q; i]) = E, we obtain 



D E = j [dq (D pp d p H + D pq d q H) - dp(D qq d q H + D pq d p H)] 

B) [ [ dqdp[d p {D pp d p H + D pq d q H)+d q (D qq d q H + D pq d p H)] (18) 
J Ja(e) 



2tt 



Near the bifurcation point the momentum p is small, p oc 7y 3 / 4 and p oc 77 1 / 2 for the local and nonlocal bifurcation, 
respectively. The derivatives d q H,d 2 H are also small, with d 2 H oc -q 1 ! 2 and d 2 H oc rj for the local and nonlocal 
bifurcation. Then taking into account that d p H = p, we can simplify Eq. H18[) . 

f[ dqdpD ppd 2 p H=^l 
IA(E) Z7T J JA(E) 



D E « / / dqdpD pp d;H = -A-± / / dqd P D pp . (19) 



Here, we have assumed that the coefficient D pp is not 
proportional to a power of a small parameter 77 inside 
the areas enclosed by the homoclinic orbits in Fig. ^ 

VI. ESCAPE RATE 

Diffusion over energy leads to escape from a metastablc 
state. The escape rate can be calculated from the one- 
dimensional Fokker-Planck equation for the energy dis- 
tribution p(E), 

dtp = d E (-V EP + D E d E p) (20) 

(here, again, we have disregarded the corrections oc D to 
the energy drift velocity). 

The escape rate W is determined by the probability 
density to have the energy of the saddle point (the top 
of the potential barrier) E = E s , given that, for small 
noise intensity, the system is mostly fluctuating about 
the center, E — E c . From Eq. (|2U)l we have 

W = const x exp(— R/D), 

rE, D 

R=- dEV E — (21) 
Je c D e 

(here, D is the characteristic noise intensity, see below). 
Eq. (|21J) is central for the following calculation, ft im- 
mediately allows us to find the scaling of the activation 
energy near a bifurcation point for the local and nonlocal 
bifurcation. 

In the case of the local bifurcation the ratio Ve/De 
can be obtained in the explicit form. From Eq. (|14f) we 
have Vv = a p + a q . Because the area of integration A(E) 



in Eq. (|19f) is limited to a small vicinity of the point q = 
p = 0, we can set in Eq. ijT§|) D pp (p, q) w D pp (0, 0) = D. 
Then from Eqs. UBJ), (UHJ) we obtain V E /D E = -e(a p + 
a q )/D. The activation energy of escape i?i oc is simply 
proportional to the energy difference AE = E s — E c Ijllfl. 

R\oc = e{a p + a q )AE = -s(a p + a q ) if' 2 . (22) 

Eq. IL'l'l) shows that, for the local bifurcation, in the 
general case of a noncquilibrium underdamped system 
the activation energy of escape scales with the distance 
to the bifurcation point as rp with £ = 3/2. This is the 
same exponent as in the overdamped regime. 

In the case of the nonlocal Hamiltonian bifurcation the 
range of integration over p in the expressions (|16fl , I19|) 
for Ve and De is narrow, oc 77 1 / 2 , whereas the range 
of integration over q is independent of rj. Therefore the 
integrands Vv and D pp in these expressions can be cal- 
culated for p = 0. Then integration over p can be done 
directly, and the integrals in Eqs. I|16f) and (|19fl take a 
form 

(E s -E c f/ 2 dqf{q) 

x^x-^- w jU(q)-U(q c )]j , (23) 

where x = x{E) = (E — E C )/(E S — E c ) is the reduced 
energy, and f(q) is proportional to Vv(p = 0, q) and 
D pp (p = Q,q), respectively. The functions q m i n and g m ax 
determine the limits of integration and are given by the 
equation E = rjU (q). 

Since AE = E s - E c = const x 77 [cf . Eq. lfT5)l] , the 
expression has the form (E s — E c ) 1 / 2 f(x), i.e., it 
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depends on energy only in terms of the parameter x(E) 
and is proportional to t] 1 ! 2 . Therefore the ratio Ve/De 
is a function of (E — E c )/ (E s — E c ), which is proportional 
to s/D, with D being the characteristic noise intensity. 
Then, from Eq. (|21() , to leading order in r\ the activation 
energy near the nonlocal bifurcation R n \ is proportional 
to the energy difference AE oc 77, as in the local case, 

R n \ = const x £7]. (24) 

From Eq. (|24J) the activation energy scales with the 
distance to the bifurcation point as 77^ with £ = 1. We 
are not aware of a prediction of such dependence in the 
general case. It is much slower than for the local Hamil- 
tonian bifurcation. 



VII. EXAMPLE: A RESONANTLY DRIVEN 
NONLINEAR OSCILLATOR 

A system where the both types of bifurcations may oc- 
cur and which has recently attracted significant interest 
0- Q- S 01 is an underdamped nonlinear oscillator driven 
by a resonant field. The motion of such an oscillator in 
the rotating frame is described by Eqs. with 

H(p,q) = \(<?+p 2 -lf-^\ 

v {q) = q, v {p) = p. (25) 

For several microscopic and phenomcnological models of 
noise, which include those of interest for the experiment 
0,Q,SI3 > t ne components f( q \ are independent and 
(5-correlated, with coordinate-independent D pp = D qq = 
D, D pq = 0. 

The parameter (3 characterizes the intensity of the driv- 
ing resonant field. For a small friction coefficient e, the 
oscillator is bistable in the region f3 { B ] < (3 < (3 { g\ where 

« e 2 and (3 ( B ] « 4/27. 

The Hamiltonian H is shown in Fig. [2 It has the shape 
of a tilted Mexican hat. The local maximum (the top of 
the internal dome) A\ and the absolute minimum Ai are 
centers; they become attractors due to friction; S is the 
saddle point. A general expression for the escape rate of 
the oscillator was obtained in Ref. Il9t in particular, the 
activation energy was calculated explicitly as a function 
of (3 in the underdamped case e<l. 

For the bifurcation value of the control parameter 

(2) 

(3 = /3 B = 4/27 in the limit of small damping, the sad- 
dle point S in Fig. |2 merges with the dome top A\ at 

= -1/V3,p = 0. For 77 = (4/27) 1 / 2 - (3 1 / 2 < 1 and 
(2) 

small Sq = q — q B , p, to the leading order in Sq, p, 77 the 
Hamiltonian H25|) has the form 

H w -^p 2 - ^(^) 3 + v(Sq) + const. (26) 

This expression has essentially the same form as the 
Hamiltonian © near the local bifurcation point. Note, 



H 




FIG. 2: The effective Hamiltonian H (p, q) of a resonantly 
driven nonlinear oscillator in the rotating frame 11251 . The 
plot refers to the reduced field (3 = 2/27. The local maximum 
Ai and the minimum A2 are the centers, and S is the saddle 
point. The trajectories for constant energy H(p,q) = E are 
shown in the lower panel. 

however, that the center is the local maximum rather 
than the minimum of the potential. 

It follows from Eq. I|26|l that near the local bifurca- 
tion point \AE\ = \E S - E c \ = (4/3 5 / 4 )ry 3 / 2 . Then from 
Eq. G3 

R loc EE Rf = (4/3 1 /4) £r? 3/2 = | (/? (2) _ ^3/2 (2?) 

(we have taken into account that, for the Hamiltonian 
()26J) . the last term in Eq. Q19|l has an extra factor 3/2). 

Eq. Q27f) coincides with the expression obtained for the 
present system by Dmitriev and Dyakonov |17| . It coin- 
cides also with the expression obtained from the analysis 

of the oscillator dynamics very close to the bifurcation 

(2) 

point p B , where the oscillator motion is overdamped 
[lOj . This is similar to the situation in equilibrium sys- 
tems, where the escape rate as a function of the distance 
to the bifurcation point is given by the same expression 
in the opposite limits of underdamped and overdamped 
motion. 

The dynamics of the oscillator near the bifurcation 
point f3 B ^ is described by the nonlocal bifurcation the- 
ory. The structure of the phase portrait in the oscillator 
variables (q,p) for the appropriate energy range is clear 
from Fig. [21 The homoclinic trajectories form a double 
loop with narrow space between the trajectories. The 
trajectories that surround the center Ai have the shape 
of horseshoes "squeezed" into the interloop space, as seen 
also from Fig. 

The slowing down of motion is a result of the fact that 
the system moves in the opposite directions along the 
homoclinic trajectories. When they approach each other 
(and the circle p 2 + q 2 = 1 that lies between them) with (3 

approaching (3 B ^ , the motion is slowed down everywhere 
between them, including the center. 
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1.5 




FIG. 3: Homoclinic trajectories of a resonantly driven oscil- 
lator near a nonlocal Hamiltonian bifurcation. The center A2 
is marked by the star,,? is the saddle point. As the control 
parameter f3 approaches the bifurcation value /3^' the homo- 
clinic trajectories approach each other and ultimately merge 
at the bifurcation point. At this time the energies of the sad- 
dle and the center become equal to each other. The data refer 
to p = 3.9 x 10" 4 . 



Even though the phase portrait in Fig. looks some- 
what different from the one in the lower panel of Fig. ^ 
they can be mapped onto each other. Therefore the re- 
sults for escape near a nonlocal bifurcation should apply 
to the oscillator. 

For j3 <C 1 and \p 2 + q 2 — 1| <C 1 it is convenient 
to change in Eq. (|25[) to new canonical variables, P = 
— (p 2 + q 2 — 1)/2 and Q = arctan(p/q). In these variables, 
to leading order in P, f3 the Hamiltonian takes the form 

©, 

H = P 2 -77 cos Q, r? = /3 1/2 . (28) 

The center and the saddle point are Q c = 0, P c = 
and Q s = 7r, P s = 0, respectively. The energy difference 
AE = E s - E c = 277. 

From Eq. I|25[) and the expressions for the noise dif- 
fusion coefficients, the corresponding functions that de- 
scribe friction force and noise in variables Q,P are 

v {p) = 1-2P, v iQ} = Dp Q = 0, 
D PP = D(l - 2P), D QQ = D(1-2P)- 1 . (29) 

We note that the oscillator has only one saddle-node bi- 
furcation point for (3 = ?y 2 <C 1, because by construction 
?? > 0. 

From Eqs. flU), (JTSJl, (|23> we obtain V B /D B = e/D. 
Therefore the activation energy of escape near the non- 
local bifurcation is 

Rni = R ( b = 2er] = 2eP 1/2 . (30) 

Eq. (|30[) coincides with the expression in Ref . obtained 
by a completely different method. 



VIII. DISCUSSION OF RESULTS 

We have considered activated escape in systems close 
to a saddle- node bifurcation point, the problem of inter- 
est not only for recent experiments in different fields of 
physics 0,01 HE EJi but also for such areas of recent at- 
tention as epidemic control ■ Of special interest to us 
were undcrdamped systems. We studied the range of the 
control parameter where a system is close, but not too 
close to the bifurcation point. In this range the motion 
is already slowed down and is described by one degree 
of freedom, but is still underdamped. The slowness of 
motion justifies describing it by a Langevin equation in 
which delay of the friction force and finite correlation 
time of noise are disregarded. 

We have identified two scenarios which an under- 
damped system can follow as it approaches a saddle-node 
bifurcation. One is local, where in the neglect of damp- 
ing the system has a saddle point and a center close to 
each other in phase space. As the control parameter 77 
approaches the bifurcation value r/s — these states ap- 
proach each other. It is the closeness of the states that 
leads to motion slowing down in this case. The second 
scenario is nonlocal, where the saddle point and the cen- 
ter remain far from each other until damping becomes 
strong. However, the phase trajectories have a shape of 
narrow cigars, and since on the opposite sides of the cigar 
the system moves in opposite directions, the overall mo- 
tion becomes slow. As r\ — > j]b, the cigars are further 
squeezed, leading to further slowing down. 

The major effect of friction is energy drift towards a 
stable state. In contrast, noise leads to energy diffu- 
sion away from the metastable state and ultimately to 
escape. The activation energy of escape R is determined 
by the ratio of the drift and diffusion coefficients. We 
found that, for both bifurcation scenarios, this ratio has 
a simple form. It can be found explicitly for the local 
scenario, whereas for the nonlocal scenario it has a form 
of a function of reduced energy. As a result, the activa- 
tion energy is proportional to the energy difference in the 
saddle point and the center. 

We have shown that, for both types of underdamped 
systems, the activation energy displays scaling R cx 77^. 
The scaling exponent £ = 3/2 for the local scenario and 
£ = 1 for the nonlocal one. The scaling £ = 3/2 has 
been known for a saddle-node bifurcation in overdamped 
systems and also for escape of a Brownian particle due 
to thermal fluctuations. As we show, this behavior is in- 
dependent of the nature of friction and fluctuations and 
occurs even in the absence of detailed balance. The scal- 
ing £ = 1 is a new result. We give an example of a 
system of recent interest where both types of scaling can 
be observed. 
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